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ABSTRACT 

We have determined the CMB temperature, T(z), at redshifts in the range 
0.023-0.546, from multi-frequency measurements of the S-Z effect towards 13 
clusters. We extract the parameter a in the redshift scaling T(z) = Tq(1 + 
z) 1 ^ , which contrasts the prediction of the standard model (a = 0) with that 
in non-adiabatic evolution conjectured in some alternative cosmological models. 
The statistical analysis is based on two main approaches: using ratios of the 
S-Z intensity change, AI, thus taking advantage of the weak dependence of the 
ratios on IC gas properties, and using directly the AI measurements. In the 
former method dependence on the Thomson optical depth and gas temperature 
is only second order in these quantities. In the second method we marginalize 
over these quantities which appear to first order in the intensity change. The 
marginalization itself is done in two ways - by direct integrations, and by a 
Monte Carlo Markov Chain approach. Employing these different methods we 
obtain two sets of results that are consistent with a = 0, in agreement with the 
prediction of the standard model. 
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Introduction 



The variation of the Cosmic Microwave Background (CMB) temperature with redshift 
is a basic relation which in adiabatically evolving cosmological models is T(z) = Tq(1 + z), 
normalized to the CO BE/FIRAS value at the present epoch, T = 2.725 ±0.002K (author?) 
(IMather et al. 19991 ). In light of the fundamental role the CMB plays in cosmology, and given 
our detailed knowledge of its spatial structure, the lack of precise observational confirmation 
of this important relation is rather surprising. Further motivation to measure T{z) is pro- 
vided by the prediction of different non-linear redshift scaling law in alternative cosm ological 
models (e.g. (author?) lOverduin fc Cooperstock 19981. iMatviasek 19951 . iPuy 2004 ). In par- 
ticular, we consider here the scaling law proposed by (ILima et al. 20001 ) . T(z) = T (l + zy~ a , 
with a a free (constant) parameter; this scaling is presumed to be a consequence of photon 
number and radiation entropy non-conservation. 

The T(z) dependence has traditionally been tested by measurements of quasar ab- 
sorption spectra whi ch include fine-structure lines from interstellar atoms and ions (e.g., 
(ILoSecco et al. 20011 )). Using this method the CMB temperature was measured to a max- 
imum redshift of 3.025, at which T ~ 12.6±^ K was determined from an analysis of 
the CH fine-structu re lines in the damped Ly a system towards the quasar Q0347-3819 
(IMolaro et al. 20021 ) . The fact that the CMB is not the only radiation field populating 
the energy levels (from which the transitions occur), and lack of detailed knowledge of 
the physical conditions i n the absorbing clouds - i s a major source of systematic uncer- 
tainty (e.g., (author? ) ICombes and Wiklind 19991 ) . A strong limit on a was deduced 



(jOpher &: Pelison 20051 ) using CMB and galaxy distribution data. However, this method 
cannot be used to determine the T(z) scaling. 

The possibility of determining T(z) from measurements of the Sunyaev-Zeldovich (S-Z) 
effect had been suggested long ago (IFabbri. Melchiorri fc Natale 19781 ; iRephaeli 19801 ) . The 
effect - Compton scattering of the CMB by hot intracluster (IC) gas - is a small change 
of the CMB spectral intensity, Al, which depends on the integrated IC gas pressure along 
the line of sight (los) to the cluster. The steep frequency dependence of the change in 
the CMB spectral intensity, Al, due to the S-Z effect allows the CMB temperature to be 
estimated at the redshift of the cluster. Since the ratio of the values of Al measured at two 
frequencies is essentially independent of the cluster properties, t he value of the temperature 
at the cluster redshift can be deduced directly from this ratio (IRephaeli 19801 ) . The much 
improved capability of precise multi-frequency measurements of the effect enhances interest 
in this method to measure T(z) in nearby and moderately distant clusters. 

The first attempt to determine T(z) from an analysis of multif requency S-Z measure - 
ments in the Coma and A2163 clusters of galaxies was reported in (jBattistelli et al. 20021 ). 
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Preliminary resul ts from an analysis of a larger sample of clusters, with z = 0.023 — 0.546, 
was presented in ( IDe Gregori et al. 20081 ). In this paper we report the final results from the 
current cluster sample carrying out a significantly expanded statistical data analysis. The 
statistical analysis is based on both ratios of AI, which are weakly dependent on the cluster 
properties, as well as the values of the individual AI whose dependence on cluster properties 
is marginalized over (also) the Thomson optical depth. We discuss the consistency between 
these two analysis methods. 

In Section [2] we outline the basic methods used to determine T(z). Data analysis and 
the results for a are presented in Section [31 We end with a brief summary in Section HJ 



2. Method 



Compton scattering of the CMB in a cluster causes a change of intensity that can be 
written as 



AI 



2k 3 T 3 x 4 e x 



dr 



eh{x)-(3 + R{x,e,(3) 



h 2 c 2 ( e *-l)2_ 

where x = hv/kT is the dimensionless frequency, 9 = kT e /mc 2 is the electron temperature 
in units of the electron rest energy, and f3 is the line of sight (los) component of the cluster 
velocity (divided b y c) in the CMB frame. The integral is over t he Thomson optical depth, r . 
Both the thermal (jSunyaev fc Zel'dovich 1972T ) and kinematic (jSunyaev fc Zel'dovich 19801 ) 
components of the effect are included in Equation (1) in the first two terms, and in the 
function R(x,0,/3). This latter term is the relativistic correction (IRephaeli 1995T ) to the 
non-relativistic expressions for the thermal and kinematic components. An analytic approx- 
imation (which is sufficiently exact even close to the crossover frequency) can be written in 
the form 

R(x,6,P) ~^[/ 2 (i) + e/ 3 (i)+eV4W+e^5W]-^[9iW+^(a:)]+/? 2 [l + ^W] . (2) 



The spectral functions fi and gt are specified in ( jShimon and Rephaeli 20041 ); see also (jltoh et al. 20021 ) 



With T(z) = T(0)(1 + *) (1-a) , and v = u Q (l + z), the non-dimensional frequency 
depends on redshift if a ^ 0, x = Xo(l + z) a , with xo = hiso/kT(0). The steep frequency 
dependence of the change in the CMB spectral intensity allows the CMB temperature to 
be estimated at the redshift of the cluster. In the non-relativistic limit AI depends linearly 
on the Comptonization parameter, y = f 9dr, which includes all dependence on IC gas 
properties. A ratio of values of AI at two frequencies is then essentially independent of the 
cluster properties. In the more general case, the first term in the square parentheses in Eq. 
([!]) still dominates over the other two, except near the crossover frequency, where the sum 
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of the temperature dependent terms vanishes. For values of x outside some range (roughly, 
3.5 < x < 4.5), the dependence of AI on j3 is very weak since the observed temperature 
range in clusters corresponds to 0.006 < 6 < 0.03, whereas typically (3 < 0.002. Therefore, 
data points at or near the crossover frequency are noise-dominated, and thus carry relatively 
low statistical weight in the likelihood function we construct for a or T(z). Nevertheless, we 
include these points in our analysis. 



3. Data Analysis 

We analyzed results of multi-frequency S-Z measurements toward 13 clusters spanning 
the redshift range 0.023 — 0.550. The dataset includes measurements with the BIMA, OVRO, 
SUZIE-II, SCUBA, and MITO telescopes, as well as gas temperatures determined from X- 
ray measurements. We assumed a gaussian profile for the spectral bands of each experiment. 
In order to evaluate systematics induced by different spectral efficiencies of the bands, we 
repeated the analysis with squared profiles. The error contribution is negligible as compared 
with uncertainties in values of the central frequency and bandwidth. The clusters, redshifts, 
and gas temperatures are listed in Table (U and the S-Z data are in Table [2j In the analysis 
we allowed for calibration uncertainty, considered as an uncertain scale factor, which was 
modeled as a gaussian with mean 1 and 10% standard deviation []. This is an adequate 
approximation given the relatively narrow widths of the respective gaussian distributions. 
For the radial component of the cluster peculiar velocities we assume a universal vanishing 
mean with a 1000 km/s standard deviation. 

The calculation begins with convolution of AJ(i/) with the detector Gaussian spectral 
response function with width a UQ , 

i poo __{"- v o) 2 

AI(v ) = —= — / AI(u)e 2CT "o dp, (3) 



2ixa UQ j u =o 

using the approximate analy tic expression for (the relativist ically correct) AI{v) derived by 



(IShimon and Rephaeli 20041 ) . Integration by parts then yielded an analytic approximation 
to the frequency-convolved AI(v), accurate to order (cx^ /z/) 4 . The sole purpose of deriving 
the latter analytic (rather cumbersome) approximations (which are not specified here) was to 
enable a faster convolution of the theoretical expression for AI(u) with the detector Gaussian 
spectral response. 



lr The calibration errors for the experiments considered are at a maximum level of 10%, based on obser- 
vations of the brightness of planets. We have adopted a conservative approach in assuming maximum errors 
in order not to be biased by different calibration procedures of different experiments. 
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The probability of measuring an intensity change AI b s (uj) at frequency Uj towards a 
cluster at redshift z i: with velocity parameter /3 h gas parameters Tj, 9i, and model parameter 
a is, 

P[r it U A; «, a) oc exp{- £ [^ft.^^ -A^M]^ (4) 

where AJ is the result of convolution of the theoretically predicted intensity change with the 
spectral response function (calculated using Eqj3]), and <Ti(vj) is measurement error standard 
deviation. 

The CMB temperature at the redshift of each cluster was first extracted by performing 
statistical analysis on the ratios = AT(z/,)/AT(z/j) (AT(z/j) is the observed thermody- 
namic temperatures in the i-th photometric band) marginalizing over 6 and (3 - the original 
(hereafter) RI approach - and by using each of these separately and marginalizing also over 
r, 6 and (3 - the DI approach 0. Assuming the S-Z data are gaussianly distributed allows us to 
construct the relevant single-cluster and multiple-cluster likelihood functions. We note that 
the high-dimensionality of the marginalization, over the full parameter space, including also 
the calibration uncertainty (see Section 13. ip can potentially be a source for large numerical 
errors. We therefore employed two separate pipelines in the DI analysis - a direct numerical 
integration and Monte Carlo Markov Chain (MCMC) sampling. After proper treatment of 
the systematics the different methods give consistent results. 



3.1. Likelihood of Intensity Ratios 



The RI approach - first employed in our previous work (jBattisteUi et al. 20021 ) - is based 
on the use of the likelihood function of intensi ty ratios, there by removing (to first order) de- 
pendence on the Comptonization parameter (jRephaeli 1980 ). By doing so we largely avoid 



the need to account for modeling uncertainties in the determination of the gas density and 
temperature profiles (from X-ray surface brightness and spectral measurements). Depen- 
dence on the measured gas temperatures, T e , is also largely removed in the leading term, 
although T e is still needed to calculate the relativistic correction which includes second order 
terms in T e . Measurements at the crossover frequency were not included in the analysis due 
to the inherent problematics associated with inclusion of very weak signals in the likelihood 
function (e.g., 'Cauchy tail', bimodality, etc), which are expected to be dominated by instru- 



2 Note that for marginalizing over a range of gas temperatures we have determined this range by using a 
beta profile for the gas density. The modeling uncertainty due to the use of this profile is negligible. 
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mental noise and systematic uncertainties, so their relative weight in the likelihood function 
would be negligible. 

The distribution of ratios may have complicated features, but if carefully applied can be 
a very powerful tool due to the reduced number of essential parameters in the problem. In our 
case the distribution of ratios of temperature changes is manifestly non-gaussian, but we are 
interested in constraining the CMB temperature (or the parameter a) rather than AT itself. 
The distribution is expected to be normal when the cluster sample is large. As shown in 
Appendix A, this approach has the unwanted feature that the expectation value of the ratio 
is biased to a degree which is of O ((AT/T) 2 ), but in the limit of very precise measurements 
the fractional error decreases as does the degree of bias. This follows from the fact that a 
uni-modal distribution function behaves as a gaussian sufficiently close to its peak, and if 
measurement errors are sufficiently small the function is sharply peaked. In other words, 
when the function is highly concentrated near the peak, any value of the parameter (for 
which the distribution function does not vanish) is in the 'gaussian- regime'. This guarantees 
that the ratio-distribution is an unbiased gaussian. However, for the small number of clusters 
in our present sample, and the low quality of some of the measurements, our results are still 
affected by this bias, but at a level smaller than the statistical uncertainty. The situation will 
improve substantially in the near future, when planned surveys with PLANCK and other 
S-Z projects will provide a very large sample of clusters. A more complete discussion of the 
ratio approach is given in Appendix A. 

The joint distribution function is derived from a set of n S-Z measurements (at n different 
frequency bands, excluding the crossover frequency) of a given cluster. In doing so we choose 
the S-Z data point with the smallest fractional error to be in the denominator of the n — 1 
ratios so as to minimize the bias (Appendix A). We construct the likelihood function of the 
simultaneously calculated n — 1 ratios for a given frequency (in the denominator) obtaining 
a single likelihood function per cluster. 

The likelihood function that encapsulates the multi-cluster information is defined as 



where the index i runs over all clusters in the sample. In the following we marginalize over 6i, 
Pi, and calibration uncertainty, denoted here by C. In the limit of sufficiently large dataset, 
the joint likelihood function is gaussian in a by virtue of the central limit theorem and in 
accord with the reasoning we specified above. The average gas temperatures and their errors, 
as obtained from X-ray observations, are summarized in Table HJ The smaller (than thermal) 
kinematic component vanishes to first order (due to motion either toward or away from the 
observer). 




(5) 
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3.2. The DI Approach 

As noted, we carried out another independent analysis based on using the measured 
values of Al themselves, rather than their ratios. This second approach was implemented 
in two separate and very different treatments. The first consisted of the construction of a 
likelihood function that incorporates the multi-cluster information, with direct numerical 2D 
integrations involved in the marginalization over r, 9, (3 and C (the integration over (3 and 
C is carried out analytically after Taylor expanding the integrand in powers of f3 and C) 

LW^UiLlTiArfi&a). (6) 

For t we use a flat (positive definite) prior; for the other parameters priors are chosen as in 
the RI approach (see Section EJ). 

The second independent treatment is based on MCMC sampling, which is very com- 
monly used in a wide class of problems, where complex dependencies between the parameters 
greatly impact the likelihood function, so that dire ct integrations could be quite tediou s and 



prone to numerical instabilities (e.g., (author?) iLewis and Bridle 20021 . IVerde 20071 ). In 
this treatment we construct a single likelihood for each cluster, in order to get directly T(z) 
at the cluster redshift, and also to check parameter degeneracies. 

The MCMC method constitutes a random sequence of realizations of the fit parame- 
ters which - when properly sampled in parameter space - tends to reproduce the posterior 
distribution for all the parameters. The technique is easily implemented into an iterative 
code that draws candidate sets of fitting parameters of chain elements from a given proposal 
distribution, and inserts the candidates into the chain if some proper statistical criterion 
i s met. The sampling a pproach we us ed is the one proposed by Metropolis and Hastings 



([Metropolis et al. 19531 ; Hastings 19701 ). According to this scheme, each candidate set of 
parameters $j is drawn from a gaussian centered over the previously accepted chain element 
i?t-i- The width of the proposal distribution is a critical parameter of the code which de- 
termines the efficiency of the sampling algorithm (i.e. the capability to sample from the 
multi-parameter posterior distribution in a given number of realizations, or acceptance rate) 
and its ability to explore the full domain of the fitting parameters. Different test runs of the 
MCMC code were needed to find the best tradeoff between acceptance rate and accuracy 
of the fit for given chain lengths, by adjusting the width of the proposal distributions. In 
the final configuration, the code was able to reconstruct the posteriors with 5% -r- 10% ac- 
ceptance rate and sufficient samplings of the parameter space were achieved in less than 10 6 
samplings. Converg ence of the MCMC runs was also tested through the Gelman-Rubin test 



( IGelman et al. 19921 ) which not only tests convergence but can also diagnose poor mixing. 



We have verified that the results are independent of the starting point in parameter space 
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by performing multiple shorter runs of the code and confirming the statistical consistency 
of the results of the single runs. Undersampling of the final chains was performed to reduce 
the self-correlation which is naturally induced in a sampled Markov Chain, and may there- 
fore affect the estimate of sample variances for the different parameters. Once we have the 
posteriors for the various parameters we can characterize the distributions in terms of the 
quantities of interest: expected value, mode, standard uncertainty, probability intervals. 

The above analysis has been performed for each cluster to derive T(z). We present the 
results in terms of T(z) in Table [3, where we give expected values and standard deviations. 
Note that T(z) at the cluster redshift is independent of the particular scaling assumed for 
the temperature (i.e. the Lima model in this case); the only assumption made is that the 
frequency obeys the standard scaling, u(z) = u (l + z) is valid. Rather, the a value is strictly 
related to the specific model under consideration. 

By studying the correlations between the variables, we checked that the main degen- 
eracies are between T(z) and r, and T(z) and (3. In particular, given the actual level of low 
precision S-Z measurements, only the correlation between T(z) and f3 is always evident. We 
verified this behaviour by simulating a dataset for a single cluster. We checked that, in order 
to reduce the impact of the degeneracy between T(z) and /3 and then to reduce the uncer- 
tainty in the determination of T(z), better knowledge of the peculiar velocity is required (see 
Fig- ©)• We used two different priors for the peculiar velocity: one gaussian with vanishing 
mean and standard deviation 1000 km/s and the other gaussian with vanishing mean and 
standard deviation 100 km/s. 

To obtain a we have performed a fit of the T(z) data points. Of course, since the 
distributions of T(z) for individual clusters are in general slightly non-gaussian and in ad- 
dition they are frequently skewed, performing a best fit as if they were gaussian introduces 
a bias in the result. Nevertheless we tested it is a good first approximation. The prior 
for a is flat in the range a G [0, 1], to account for the theoretical constraints of the model 
flLima et al 200oh . 

3.3. Results 

To obtain the results we report below we assumed the flat prior a G [0, 1]. Employing 
the RI approach we deduce the most probable value a = 0.0241q^; all errors in this section 
are at 68% confidence. As discussed earlier, the use of intensity ratios necessarily introduces 
bias in the inferred parameter, although we attempted to minimize the bias by selecting the 
measurement with the smallest fractional error to be in the denominator of the ratios for 
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each cluster. The more precise the measurements, the weaker is this bias. To estimate the 
bias in this result we repeated the calculation with observational errors in values of Al (used 
in the denominators) reduced to very small (relative) levels, thereby isolating the impact of 
this systematic bias. Doing so yields a = 0.038+Q;Qgg. Thus, when compared with the value 
deduced from the actual data, a = 0.024+q;q24, a bias level of ~ —0.014 is indeed a relatively 
small fraction of the statistical uncertainty, ~ 0.057. 



The analysis was repeated employing the DI approach, adopting a flat r prior in the 
interval [0, 0.05]. In the first direct 3D integrations treatment we obtain a = 0.026^oo26- The 
final alpha value we get by fitting the T(z) data points obtained with the MCMC treatment 
is a = 0.062^Q Qg2- m order to check the efficiency of this method we have also considered the 
T(z) as deduced from line transitions observations (mainly estimated by UV observations of 
interstellar clouds) (see Fig. [TJ. The improvement on the constraints on a is not substantial 
in spite of the wider redshift range (as compared to the S-Z data). The final alpha value 
we get by fitting the T(z) data points is a = 0.0411q!o4i- Due to the shape of the posterior 
(see Fig. |2j) as a consequence of the strong prior on a, it is more meaningful to present this 
result as just an upper limit, i.e. a < 0.079 at 68% probability level. 

Our three different treatments yield consistent results. More important, no significant 
evidence is found for a deviation from the redshift dependence of the CMB temperature 
predicted in the standard model. 

In principle, the intensity-ratio approach to determine T(z) would seem preferable over 
the DI method since dependence on r is only second order, as compared with first order 
linear dependence on this quantity, which therefore needs to be marginalized over in the DI 
method. This is so even if the former method is somewhat biased due to the arbitrariness 
in selecting the intensity change used in the denominator of the intensity ratios. Our RI 
upper limit, a < 0.092 at the 68% confidence level, is weaker than that obtained in the DI 
analysis, 0.059; this reflects the fact that the 68% confidence region is wider than the usual la 
uncertainty region when the distribution is non-gaussian, such as the one we get with the RI 
approach. With more precise measurements, or - equivalently - a larger number of clusters, 
the distribution of ratios would approach a gaussian and would therefore be less affected by 
this bias. Ongoing and near future surveys, including those planned with ACT, APEX-SZ, 
SPT, a nd Planck, as well as detailed S-Z mapping of a sa mple of nearby clu sters with the 



MAD flLamagna et al. 20021 1 and balloon-borne OLIMPO (IMasi et al. 20031 1 experiments, 
will provide much more precise and uniform datasets that will largely remove the bias in 
the RI approach. While there is no such bias in the DI approach, the marginalization over 
values of r lowers the usefulness of this approach. On the other hand, spatially resolved 
spectroscopic observations of galaxy clusters (as proposed with the SAGACE satellite, a 



small mission project approved for Phase A by the Italian Space Agency) would allow for 
breaking the degeneracy between T(z) and cluster parameters. 



4. Summary 

By adopting a more general redshift scaling law that includes the prediction of the 
standard model when a = 0, we are able to test both adiabatic and non-adiabatic models. 
While we obtain somewhat different results for a in the above three treatments, the degree 
of overlap of the respective uncertainty intervals implies that these values are statistically 
consistent. Already with the current sample of 13 clusters with medium-quality S-Z data, 
we are able to verify the standard scaling law at a good level of precision. More precise 
measurements of the T(z) scaling law could possibly have interesting ramifications on setting 
constraints on the variation of fundamental constants over cosmological time. 

We wish to thank G. D'Agostini for useful discussions and the referee for helpful sug- 
gestions. MS gratefully acknowledges useful discussions with Ran Shimon. This work has 
been supported in part by MIUR/PRIN 2006 (prot. 2006020237) and University of Rome, 
La Sapienza (Ateneo prot.C26A0647AJ). 



Appendix A: Ratio Distribution 

We briefly summarize some of the basic elements of the theory of ratio distributions 
employed in this work. If two variables x and y are drawn from distribution functions S(x) 
and T(y), respectively, then the probability that their ratio is in the interval [r, r + dr] is 

pod n[r+dr)x 

P(r < x/y < r + dr) = / / S(x)T(y)dx dy, (1) 

J x=—oo J y=rx 

which can be shown to imply that 

/oo 
\x\S(x)T(rx)dx (2) 
-oo 

where in the last step we took the absolute value of the Jacobian so that the distribution 
function is non-negative. For gaussian distributions S and T and finite r, this integral can 
be carried out analytically. 

In the case of three independent, gaussianly distributed, data points x, y, z, and con- 
sidering the two ratios r\ = y/x, and r 2 = z/x , where r\ = Ri(x,y) and r 2 = R 2 (x,z), the 



joint density function P(r\,r 2 ) is 



P(r 1 ,r 2 ) = J 5[n - y)]5[r 2 - Rq(x, z)]f(x, y, z)dxdydz 



(3) 



where f(x, y, z) = S(x)T(y)U(z) because data are independent by definition, and which once 
integrated, after employing y = r±x and z = r 2 x, gives 



P(r 1 , r 2 ) oc / x 2 e 



(r 1 x-p 2 ) 



'1 e 



dx 



(4) 



(see for example (author?) iD'Agostini 2003bl ). 



In the case of n ratios of gaussianly distributed n + 1 data-points we similarly obtain 
the probability for the nontrivial n ratios (in addition to the ratio r which equals unity) to 
be simultaneously r 1; r 2 ,....,r n 



P(r u r 2 , ...,r n ;9,p r ;a) oc 



dx ■ \x\ n exp 



n+l 

-E 



(xrj_i - Pi)' 

2a? 



(5) 



where pi and <jj are the expectation values and widths of the individual gaussians. In our 
case (Table 2) n can be either 1, 2, or 3. We carried out these ID integrations numerically 
In practice, the ratios r« are the theoretical ratios of SZ observations at the various frequency 
bands which are functions of the gas temperatures 6 and cluster peculiar velocities f3 r , as 
well as the CMB temperature or a. The values pi are measured with errors <r, (Table 2). 
The next step is to marginalize Eq. (A5) over the gas temperature 9 and cluster velocity f3 r . 

A generic problem of ratio distributions is their inherent bias, which can be illustrated 
as follows. Assume two measurements p\ and p 2 with gaussian errors a% and a 2 , respectively. 
The probability function for the ratio r is 



P(r) 



2no~iO~ 2 



xe 



T i e 



dx. 



(6) 



It can be easily verified that this probability function is indeed normalized, J™ P(r)dr = 1. 
Similarly, we can calculate the expectation of the ratio r; (r) = P(r)rdr. A change of 
variables leads to 



1 



27r<7i<72 



dx 



e 2a i 



(z-p 2 r 



zdze 



(7) 



The second integral can be readily calculated and yields a 2 p 2 y27r, and the first integral is 



12 



Assuming that the measurement at the denominator of the ratio is very accurate, i.e. that 
°\ ^ Pi) the denominator in the integral can be expanded in powers of and the integration 
is then carried out term by term. The leading order expectation value for the ratio is then 



This demonstrates the bias in ratio distributions; it is dominated by the error of the denomi- 
nator and therefore we choose in all our calculations the denominator to be the measurement 
with the smallest fractional error as to minimize the bias. In higher dimensional distribution 
functions the definition of the bias is not unique and therefore cannot be unambiguously 
calculated and corrected for. Note, also, that the bias is in the expectation value of the ratio 
of SZ temperatures, and not directly in a. The exact change in a as a result of this bias 
is difficult to predict a priori, but for high sensitivity measurements with large samples of 
clusters the effective fractional error drops significantly, and the bias in a is much weaker. 
As in the DI method our likelihood function is obviously non-gaussian in a and we define 
the average a as the median value of that distribution. 
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0.01 



0.10 



1.00 



Fig. 1. — T(z) as deduced from S-Z spectra from complete d ata set of clusters (A) together 
with results deduced from line transition s observations (0): (ICui et al. 20051 ; iGe et al. 19971 ; 
Srianand et al. 2000 ; iMolaro et al. 2002l ). The solid line is the best-fit to the Lima scaling. 
The dot-dashed line is the standard scaling. 
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a 

Fig. 2. — Posterior of the a parameter, as obtained by performing a fit of the T(z) data 
points. 

Table 1: Cluster redshift and gas temperatures 



Cluster 




z 


kT° (keV) 


A1656 





023 


8.25±0.10 b 


A2204 





152 


11.53 ±2.80 


A1689 





183 


9.59 ±2.80 


A520 





200 


8.33 ±0.76 


A2163 





202 


16.18 ±3.86 


A773 





216 


6.62 ± 1.30 


A2390 





232 


10.13 ± 1.22 


A1835 





252 


13.45 ±4.02 


A697 





282 


10.60 ± 1.13 


ZW3146 





291 


8.12 ± 1.00 


RXJ1347 





451 


13.69 ±2.64 


CL0016+16 





546 


10.10 ±2.57 


MS0451-0305 





550 


10.68 ±2.93 



a (|Bonamente et al. 20061 ) 
b |Arnaud et al. 200ll) 
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Fig. 3. — Parameter correlations for a simulated cluster: the contours show the 68% and 95% 
confidence limits from the marginalized distributions. Black contours are obtained allowing 
for a peculiar velocity prior with vanishing mean and standard deviation 1000 km/s; cyan 
contours are obtained allowing for a peculiar velocity prior with vanishing mean and standard 
deviation 100 km/s. T* MB =T(z)/(l + z) 
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Table 2: SZE measurements of 13 clusters by different experiments expressed in central 
thermodynamic temperature. 





OVRO+BIMA Q 




SuZIE II b 




SCUBA C 


Cluster 


AT 30 GHz 


AT145GHZ 


AT 2 2ighz 


AT273GHZ 


AT355GHZ 


AT353GHZ 




(mK) 


(mK) 


(mK) 


(mK) 


(mK) 


(mK) 


A520 


-0.66±0.09 d 


-0.44 ± 0.13 


0.14 ±0.14 




1.78 ± 0.41 


2.60 ±0.57 


A697 


-1.22 ±0.12 


-0.93 ±0.13 


0.41 ±0.16 




2.80 ±0.62 




A773 


-1.08 ± 0.11 


-0.91 ±0.16 


0.04 ±0.25 




2.40 ±0.89 


2.7±2.0 


A1689 


-2.06 ±0.17 










2.93 ±0.40 


A1835 


-2.90 ±0.21 


-1.74 ±0.26 


0.14 ±0.41 


1.73 ±0.59 






A2204 


-3.22 ±0.32 


-0.65 ±0.10 


0.21 ±0.10 




1.85 ±0.44 




A2390 




-0.91 ±0.10 


-0.10 ±0.17 




1.23 ±0.34 


2.40 ±0.44 


CL0016+16 


-1.44 ±0.09 


-0.57±0.23 


0.66 ±0.46 


1.82 ±0.68 




1.96 ±0.64 


MS0451-03 


-1.48 ±0.09 


-0.779 ±0.065 


-0.21 ±0.10 


0.91 ±0.66 


1.16 ±0.34 




RXJ1347 


-5.15 ±0.60 


-3.22 ±0.39 


-0.10 ±0.39 




6.2 ± 1.5 


5.36 ±0.54 


ZW3146 


-2.02 ±0.25 


-1.56 ±0.39 


-0.25 ±0.48 




3.1 ±1.3 






OVRO d 




MITO e 






SCUBA C 


Cluster 


AT32GHz 


AT143GHZ 


AT214GHZ 


AT272GHZ 




AT353GHz 




(mK) 


(mK) 


(mK) 


(mK) 




(mK) 


A1656 


-0.520 ±0.093 


-0.179 ±0.037 


0.033 ±0.080 


0.170 ±0.034 








OVRO+BIMA a 




SuZIE P 






SCUBA C 


Cluster 


AT30GHZ 


AT142GHZ 


AT 2 i7GHz 


AT 2 68GHz 




AT353GHZ 




(mK) 


(mK) 


(mK) 


(mK) 




(mK) 


A2163 


-1.89 ±0.17 


-1.011 ±0.098 


-0.21 ±0.16 


0.66 ±0.24 







Bonamente et al. 2006 


): b ( Benson et aZ. 2003: Benson et al. 2004): 


Zemcov et al. 20071): d (Herbie et al. 1995; Mason et al. 2001): 


De Petris et al. 2002; 


Savini et al. 2003)/ ( 


Holzaofel et al. 1997a) 
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Table 3: T(z) values as estimated for each 



deviations. 


Cluster 


T(z) 




(K) 



A1656 


2. 


.72 


± 





,10 


A2204 


2 


,90 


± 


0. 


.17 


A1689 


2. 


.95 


± 





,27 


A520 


2 


.74 


± 


0. 


.28 


A2163 


3 


.36 


± 





,20 


A773 


3 


.85 


± 


0. 


.64 


A2390 


3 


.51 


± 


0. 


.25 


A1835 


3 


.39 


± 


0. 


.26 


A697 


3 


.22 


± 





.26 


ZW3146 


4. 


.05 


± 


0. 


.66 


RXJ1347 


3 


.97 


± 





.19 


CL0016+16 


3 


.69 


± 





,37 


MS0451 


4 


.59 


± 


0. 


.36 



. We report expected values and standard 



